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Abstract 
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PACS: 02.30Hq, 02.301k, 02.20Sv, 45.40Cc 

Keywords: Lie symmetries, first integrals, Kowalevski top 

Running title: Lie symmetries & first integrals: Kowalevski top 



Corresponding author. E-mail: nucci@unipg.it 



1 



1 Introduction 



In January 2001, the first Whiteman prize for notable exposition on the 
history of mathematics was awarded to Thomas Hawkins by the American 
Mathematical Society. In the citation, published in the Notices of AMS 48 
416-417 (2001), one reads that Thomas Hawkins ". . . has written extensively 
on the history of Lie groups. In particular, he has traced their origins to 
[Lie's] work in the 1870s on differential equations . . . the idee fixe guiding 
Lie's work was the development of a Galois theory of differential equations 
. . . [Hawkins's book ^U]] highlights the fascinating interaction of geometry, 
analysis, mathematical physics, algebra, and topology . . .". 
In the Introduction of his book |36j . Olver wrote that "it is impossible to 
overestimate the importance of Lie's contribution to modern science and 
mathematics. Nevertheless, anyone who is already familiar with [it] ... is 
perhaps surprised to know that its original inspirational source was the field 
of differential equations" . 

Lie's monumental work on transformation groups [20], |22j . and in par- 
ticular contact transformations [23], led him to achieve his goal 24 . Lie 
group analysis is indeed the most powerful tool to find the general solution 
of ordinary differential equations. Any known integration technique can be 
shown to be a particular case of a general integration method based on the 
derivation of the continuous group of symmetries admitted by the differential 
equation, i.e. the Lie symmetry algebra. In particular, Bianchi's theorem 
0, jHH], states that if an admitted n-dimensional solvable Lie symmetry 
algebra is found, then the general solution of the corresponding n order sys- 
tem of ordinary differential equations can be obtained by quadratures. The 
admitted Lie symmetry algebra can be easily derived by a straightforward 
although lengthy procedure. As computer algebra softwares become widely 
used, the integration of systems of ordinary differential equations by means 
of Lie group analysis is getting easier to carry out. 

A major drawback of Lie's method is that it is useless when applied to sys- 
tems of n first order equations, because they admit an infinite number of 
symmetries, and there is no systematic way to find even one-dimensional 
Lie symmetry algebra, apart from trivial groups like translations in time 
admitted by autonomous systems. One may try to derive an admitted n- 
dimensional solvable Lie symmetry algebra by making an ansatz on the form 
of its generators. 

However, Nucci |30j has remarked that any system of n first order equations 
could be transformed into an equivalent system where at least one of the 
equations is of second order. Then, the admitted Lie symmetry algebra is 
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no longer infinite dimensional, and non trivial symmetries of the original 
system could be retrieved jSOJ. This idea has been successfully applied in 
several instances (SO], @3], 0, [20], |H2], EH], EH], EH- 
Here we show another striking application of such an idea. If we consider 
a system of first order equations, and by eliminating one of the dependent 
variables derive an equivalent system which has one equation of second or- 
der, then Lie group analysis applied to that equivalent system yields the 
first integral(s) of the original system which does not contain the eliminated 
dependent variable. Of course, if such first integrals exist. The procedure 
should be repeated on as many times as there are dependent variables in 
order to find all such first integrals. The first integrals correspond to the 
characteristic curves of determining equations of parabolic type |42| which 
are constructed by the method of Lie group analysis. 

We would like to remark that interactive (not automatic) programs for cal- 
culating Lie point symmetries such as [2H], [22] are most appropriate for 
performing this task. 

It is well known that if one finds a transformation which leaves invariant 
a functional describing a variational problem, then Noether's theorem [22] 
provides a first integral of the corresponding Euler-Lagrange system. Unfor- 
tunately, a general method for finding such a transformation does not exist. 
In addition, many equations of physical interest (e.g., Lorenz system in me- 
teorology j2S]) do not come from a variational problem. On the contrary, 
our method can be applied to any system of ordinary differential equations, 
even if they do not derive from a variational problem an d we do not 
make any a priori hypothesis on the form of the first integrals, apart missing 
one of the unknowns. 

In the next section, we describe the method in detail, in sections 3 and 4, 
we present the classical example of the Kowalevski top, and in section 5 the 
three-dimensional Kepler problem in cartesian coordinates. The last section 
contains some final comments. 

2 Outline of the method 

Let us consider the following autonomous (which could also be non-autonomous) 
system of N first order ordinary differential equations 




(1) 



w N = F N (wi, w 2 , • • • ,W N ) 
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Let 



/ = I(wi,w 2 , ■ ■ .,w s -i,w s+1 , . . . ,w N ), 



(2) 



be a first integral which does not depend on w s , and 

N 

X = V(t,wt, . . . ,w N )d t + G k (t,wi, . . .,w N )d Wk (3) 

fc=i 

be a generator of a Lie point symmetry group for If we derive w s from 
one of the equations say the first, then we obtain a system of N — 2 
equations of first order in W2, ■ ■ ■ , w s -i, w s+ i, . . . , wn and one of second order 
in w\. We remark that the method does not depend on the equation we 
choose from (^Q) to derive w s . After introducing the new notation Uj , (j = 
1, . . . , N — 1), we can write the system we obtain as 

U\ = /l(ui,U 2 , • • • ,UN-l,Ul) 
U2 = f2(ui,U2, ■ ■ ■ ,UN-l,Ul) 

■ (4) 
-iijv-l = fN-l(ui,U 2 ,...,UN-i,Ui) 

A generator of a Lie point symmetry group for @ is 

N-l 

X = V(t,ut,. . .,u N -i)d t + Y Gj{t,ui, . . .,u N _i)d U] (5) 

i=i 

If we apply Lie group analysis to system Q using the interactive REDUCE 
programs developed by Nucci j^H] ; , then we obtain a determining equa- 
tion of parabolic type for V. Its characteristic curves will yield m < N — 1 
transformations, which eliminate u\ from all the first order equations in 
(|H>. Thus, we have obtained a system of N — 2 equations of first order and 
one equation of second order in the new dependent variables Uj such that 
u\ = u\ and each of the other variables Uj are either the original Uj itself, 
if ii\ did not appear in the j-equation of system (@J), or the corresponding 
characteristic curve. If we apply Lie group analysis to this final system, 
then again a determining equation of parabolic type will be derived, and its 
characteristic curve, when rewritten in the original variables, will be exactly 
the first integral ©. 

Now let us consider a system of M second order ordinary differential equa- 
tions 

Xi = Hi{xi, . . . ,x M ,xi, ■ ■ ■ ,x M ), (i = l,...,M). (6) 
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A generator of a Lie point symmetry group for this system has the form 



M 



r = r(t,xi, . . .,x M )d t + ^2iji(t,xi, . . .,XM)d Xi . 



(7) 



i=i 



System (JEJ) can be converted into the following autonomous system of 2M 
first order ordinary differential equations 



At this point, we could either proceed as indicated above or choose one of 
the dependent variables to be the new independent variable y in order to 
reduce the order of system Q by one [3U]. For example, we could take 
xm = wm = V- Then, system (JHJ) becomes the following non-autonomous 
system of 2M — 1 first order ordinary differential equations with independent 
variable y 



where {h = 1, . . . , M — 1). Now, our method can be applied to this system 
as if it was system (0) . The fact that system © is not autonomous does not 
effect the result, as we will show in the case of the three dimensional Kepler 
problem in cartesian coordinates. 

The same method can be applied to a single ordinary differential equation 
of order N which can be easily transformed into a system of N equations of 
first order. It should be noticed that there could be several different ways 
of transforming an equation of order N into a system of iV equations of first 
order. Then, the just described method may give different results, videlicet 
(viz) no first integrals with certain reductions, all the first integrals with 
different reductions. 

3 Finding the Kowalevski top 

The motion of a heavy rigid point about a fixed point is one of the most 
famous problems of classical mechanics |7]. In 1750, Euler jHj derived the 
equations of motion which now bear his name, and described what is nowa- 
days known as the Euler-Poinsot case because of the geometrical description 




(8) 




— W M +h = H h (wi,...,WM-l,y,W M +l,---,W2M)/w 2 M 



(9) 
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given by Poinsot about hundred years later [HSJ. It was Jacobi ^3] who in- 
tegrated this case by using the elliptic functions which he had developed 
(along with Legendre, Abel and Gauss 26 ) and mastered 14 (we have 
translated this fundamental text into Italian and extensively commented 
[IT]). Another case was described by Lagrange \T3\, and it is known as the 
Lagrange-Poisson case, due to the extensive study done later by Poisson 
. This case can also be integrated by using Jacobi elliptic functions |44| . At 
the time, it seemed that other cases could easily be found and similarly in- 
tegrated. In 1855, the Prussian Academy of Science proposed this topic for 
a competition, but nobody applied 4 . The problem was so elusive that the 
German mathematicians called it the mathematical mermaid (die mathema- 
tische Nixe) |17j . More than thirty years elapsed before the Bordin prize was 
awarded to Kowalevski for finding and reducing to hyperelliptic quadratures 
the third case ^S] which is since known as the Kowalevski top. She solved 
the problem by looking for solutions which are single-valued meromorphic 
functions in the entire complex plane of the variable t |7|. Her method 
became what is now known as the Painleve-Kowalevski (or just Painleve) 
method [T2] . 

Hawkins had established "the nature and extent of Jacobi's influence upon 
Lie" [0] • It is a remarkable coincidence that the mathematical mermaid can 
also be found by using Lie group analysis as we show in the following. 
The Euler-Poisson equations describing the motion of a heavy rigid body 
about a fixed point are ^B] 

P = ({B -C)rq + mg((3z G -'yy G ) s )/A 
q = \AC - A)pr + mg(^x G - az G )\ / B 

< r = y(A — B)pq + mg(ay G — (3x G )J / C (iq) 
a = (3r — 7g 
$ = jp — ar 
7 = aq-f3p 

with A,B,C the principal moments of inertia , p(t), q(t),r(t) the components 
of the angular velocity, m the mass of the body, g the acceleration of gravity, 
x G ,y G ,z G the coordinates of the center of mass, and a(t), (3(t), ^{t) the 
component of the unit vertical vector. There are three first integrals for 
system (jlOj) : conservation of energy, i.e. 

h = - (Ap 2 + Bq 2 + cr 2 ) + rag [x G a + y G (3 + z Gl ) (11) 
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conservation of the vertical component of the angular momentum, i.e. 

h = Apa + Bqf3 + Cr-y (12) 
the length of the unit vertical vector, i.e. 

/ 3 = a 2 + /3 2 +7 2 (=l) (13) 

If we apply our method to system then we find only the first integral 
of the unit vertical vector which has p, q, r as missing variables. Kowalevski 
found that if one imposes the following conditions on the parameters: 

(1) A=B=2C 

(2) z G = 0, and either x G 7^ or y G 7^ 

then there exists a fourth integral, i.e. 

' 2 2 x G a-y G (3\ 2 ( x G (5 + y G a\ 2 
p -q -mg - J + [2pq - mg - I (14) 

We notice that 7 and r are missing in IJ14J) . Thanks to our method, we 
can find the Kowalevski top by searching for a first integral which does not 
contain 7. First we derive 7 from the second equation of system (|l()jl. i.e. 

Bq + (A — C)pr + mgz G a 

7 = 

mgx G 

which implies that x G must be different from zero. We obtain the following 
system of four equations of first order, and one of second order: 

ui = u\ {Amz G + (A - C)u 3 y G ) /Ax G 

-(A-B)(A- C)u x ul/BC + (A - C) 2 y G u 2 u 2 3 /ABx G 
+ (A - C)u lU 2u 3 z G /Bx G -(A- C){B - C) Ul ul/AB 
-(Au 2 x G - Cu 3 z G )(A - C)mgy G u A /ABCx G 
+ (A(A - 2C)u 2 x G + C(C - 2A)u 3 z G ) mgu 5 /ABC 
+{x G + z G )mgu\u^j Bx G 



(15) 



u 2 = -u 1 By G /Ax G + U2,{{C-A)y G u 2 + {B-C)x G u 1 )/Ax G 
+mgz G (-U4y G + u 5 x G )/Ax G 



(16) 



n 3 = ((^4 - B)uiu 2 + mg(u 4 y G - u 5 x G )) /C (17) 

(18) 



U4 = —iiiBui/mgx G + {C — A)uiu 2 uz/mgx G 
+(u 3 u 5 x G - uiu A z G )/x G 
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"5 

with 



u\Bu2/mgxG + {A — C)u2U^/mgx G + z G U2U^j x G — U3U4 (19) 



U\ =q, u 2 = p, u 3 = r, u 4 = a, u 5 = (3 (20) 
Now we apply Lie group analysis to system H15j) - (|19j) . An operator T 

5 

T = V(t,ui,U2,u 3 ,U4,u 5 )d t + G k (t,ux,u 2 ,Uz,u±,U5)d Uk (21) 

k=l 

is said to generate a Lie point symmetry group if its second prolongation 



r 

2 



k=l 



dt 



dV 



dt 



Uk-r- du, + — 



d /dGi 



dV\ 



dV\ 



dt V dt 



"1— ) - fi i d^J 5 " 1 



dt 7 



applied to system (|15) ) -t(19 |> . on their solutions, is identically equal to zero, 
i.e. 



USD 
CS3) 

m 
m 



n3-m 

CG3-CEJ 

n3-m 



(22) 



n3-m 

The five determining equations (1221) constitute an overdetermined system 
of linear partial differential equations in the unknowns V,Gk(k = 1,5) In 
fact, they are polynomials in ui, each coefficient of which must become 
identically equal to zero. In particular, the first determining equation in 
l|22[) is a polynomial of degree three for u\. The coefficient of highest degree 
yields an equation of parabolic type for V in four independent variables 
Ui,U2,U4,u 5 , i.e. 



A'Z 2 2 

A mg x G 



d 2 v 



2ABm 2 g 2 x G y G 



d 2 v 



du\U2 

d 2 V d 2 V 
+2A 2 Bmgu2X G - h B 2 m 2 g 2 y G -^- i; + 2AB 2 mgu 1 y G 



2A 2 Bmgu\x G 



0U1U5 



d V d V 

-2AB 2 mgu 2 y Gir — + A 2 B 2 u\-^- 2 - - 2A 2 B 2 uiu 2 



8 2 V 
duiUi 

d 2 v 

3U2U\ 

d 2 v 



(23) 



+A 2 B 2 u 2 2 



' du 2 u 5 

d 2 v 

du 2 



duj 
dV 



0M4M5 



dV 

A^Bmgxcj- AB 2 mgy G - — 

0U4 0U5 
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Its three characteristic curves yield the following transformations 

Buxyc . Bqy G 
u 2 = s 2 -. viz p = s 2 

Ax G axq^ 



Bui Bq 2 
«4 = S4 viz a = Si 



2mgxG 2mgx G 
By G m + 2Ax G u 2 . t-j By G q + 2Ar G p 

u 5 = s 5 + Bm — g V1Z P = S5+Bq — g 

2Amgx G AAmgxQ 

(24) 

with S2; S4, and S5 new unknown functions of t. 

As outlined in section 2, transformations (|24j) eliminate ui from all the first 
order equations in system ()15 [) -(|19 [) . 
In fact, system (|15 j) -(jl9 j) becomes: 

Hi = (-6A 2 BCuiu 2 x G 2 - 4:A 2 BCuix G u 2 u 3 z G 

—2A 2 Bmguix G 2 y G U5 — 2A 2 C 2 u 3 x G u\u 2 z G 
+2A 2 Bmguix G y G 2 Ui - 3A 2 'BCu x 2 y G u 3 z G 
+2A 2 BCuiu 3 x G y G - 5A 2 BCui 2 y G u 2 x G + 3AB 3 u 1 3 y G 2 
-2A 2 BCu x u 3 2 y G 2 - A 2 BC Ul 3 x G 2 - 2A 2 BCu x u 3 2 x G 2 
-2AC 3 u 3 2 x G 2 u\ + 2A 2 Cmgx G uiu A z G 2 + 2A 2 Cmgx G 3 uiu 4: 
—4A 2 Cmgx G 2 u 3 z G U5 — 4A 2 Cmgx G 3 u 2 u^ 
+2A 2 Cmgx G u 3 z G U4_y G - 7A 2 B 2 y G u l 2 u 2 x G 
+2A 2 Cmgx G 2 u 2 U4y G — AA 2 C 2 u 3 2 x G u 2 y G 
-3A 2 B 2 y G 2 Ul 3 - B 2 C 2 y GUl 2 u 3 z G - 2 Ul A 4 u 2 x G 2 
+2A 3 u 2 x G Cu 3 2 y G + 2A 3 u 2 x G Cu 1 u 3 z G + 2A 3 u\x G 2 Cu\ 
-2 Ul BC 3 y G 2 u 3 2 + 5A 3 u 2 x G B Ul 2 y G - 2A 3 u 2 x G 2 mgu 4 y G 
+4A 3 u 2 x G 2 Bui + 4ABCmguix G 2 y G U5 
—2ABCmgu\x G y G 2 u\ — 2ABC 2 ii\u 3 x G y G 
+2ABC 2 u lXG u 2 u 3 z G + 3ABC 2 Ul 2 y G u 3 z G 
+AABC 2 u lU3 2 y G 2 + 2ABC 2 uiu 3 2 x G 2 + 2AB 2 Cu x 2 y G u 3 z G 
+WAB 2 Cu 1 2 y G u 2 x G + 2AC 3 u 3 2 x G u 2 y G 
+2AC 2 mgu 3 x G 2 z G u 5 - 4B 3 Cy G 2 ui 3 
-2AC 2 mgu 3 x G z G u 4 y G + 2A 2 C 2 u 3 2 x G 2 u 1 
-A 2 BCu x 3 z G 2 + 3AB 2 C Ul 3 y G 2 
+2A 2 BCu 1 x G u 1 z G + 2A 3 u 2 x G 3 mgu 5 ) /2A 2 BCx G 2 

(25) 



u 2 = {-2A 2 u 2 u 3 x G y G + ABui 2 y G z G + 2ABu 1 u 2 x G z G 
+2ABu 1 u 3 x G 2 + 2ABu 1 u 3 y G 2 - 2ACu\u 3 x G 2 
+2ACu 2 u 3 x G y G - 2AmguiX G y G z G + 2Amgu 5 x G 2 z G 
-2BCu lU3 y G 2 )/2A 2 x G 2 
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u 5 



(2A 2 u\U2Xg — 3ABui 2 i/g — ^Bui^Axq + ^Amgu^XGUG 
-2u 5 x G 2 Amg + 3B 2 y G u 1 2 )/2ACx G 

(-2A 2 uiu 2 U3Xg + ABu^zg + 2ABui 2 U3Ug 
+2AB111U2U3XG + 2AC111U2U3XG — 2Amgu\UiXGZG 
+2Amgu 3 u 5 x G 2 - B 2 Ul 2 u 3 y G - 2BCu 1 2 u 3 y G ) /2Amgx G 2 

(2A 2 u 2 u 3 xg — ABu\ 2 U2Zg + ABu\ 2 u 3 xg — 2ABu\U2U 3 yG 
—2ACU2U3XG + 2Amgu2U4XGZG — 2Amgu 3 u/±XG 2 
-2B 2 ui 2 u 2 z G - 2B 2 u 1 2 u 3 x G + 2BCu\ 2 u 3 x G 
+2BCuiU2U 3 y G - 2Bmgu 1 u 5 x G z G ) /2Amgx G 2 



with 



U2 = 32, n 4 = S 4 , U 5 

We now apply Lie group analysis to system (J2 
f 



S5 



29]) . An operator T 



V(t,ui,U2,u 3 ,U4,u 5 )dt + Gi(t,ui,u 2 ,u 3 ,U4,u 5 )d Ul 
+G 2 (t, ui,u 2 ,u 3 ,u 4 , u 5 )du 2 + G 3 (t, ui,u 2 ,u 3 ,u 4 ,u 5 )d U3 
+Gi(t, ui,u 2 ,u 3 ,u 4 , u 5 )du 4 + G 5 (t, ui,u 2 ,u 3 ,u 4 , u 5 )du 5 



(27) 



(28) 



(29) 



(30) 



(31) 



is said to generate a Lie point symmetry group if its second prolongation 
applied to system (|25| ) -(|29 |) . on their solutions, is identically equal to zero, 



i.e. 



(E3) 



123-121 
123-121 
123-121 
123-121 



(32) 







123-121 

The five determining equations (|32j) constitute an overdetermined system 
of linear partial differential equations in the unknowns V,Gk(k = 1,5) In 
fact, they are polynomials in ui, each coefficient of which must become 
identically equal to zero. In particular, the fifth determining equation in 
(|32jl is a polynomial of degree one for u\. We call its two coefficients c5kl 

dV 

and c5k0. For the sake of simplicity, we assume Gk = 0, -r— = 0. Then, the 

at 
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coefficient of degree one, i.e. c5kl, yields 

dui 

Now, c5/c0 is a polynomial of degree five in u\. Therefore, its coefficients, call 
them c5m5, c5m4, c5m3, c5m2, c5ml, c5m0, must become identically equal 
to zero. The coefficient of degree five in u±, i.e. c5m5, yields 

—A 2 B 2 Cz G {-(A + 2B)u 2 z G + {A-2B + 2C)u 3 x G ) = (33) 
0U4 

which gives the condition on the parameter 

z G = (34) 

Then, the coefficient of degree four in u\, i.e. c5m4, yields 

/ dV dV 

-h(A- B)—mgy G x G - (A - 2B + 2C)—Cu 3 x G 

-(2A -B- 2C)—Cu 3 y G ) (A - 2B + 2C)AB 2 u 3 x G = (35) 

which gives the condition on the parameters 

A = 2B - 2C (36) 
Then, the coefficient of degree three in u\, i.e. c5m3, becomes 

— mgx G -— Cu 3 \ (2B-3C)(B-C)(B-2C)B 2 u 2 u 3 y 2 G = (37) 

which gives the further condition on the parameters 

B = 2C (38) 

Thus, we have found the Kowalevski top. We also notice that either condi- 
tion 2B = 3C or B = C leads to the Lagrange top. Finally, we are left with 
two linear first order partial differential equations in V = V(u2,u 3 ,U4,u^), 
the coefficient of degree two in m, i.e. c5m2, 

( dV dV \ dV o , o 2 , dV 

2 [d^ XG ~ d^ VG Cu ^-^ m 9x G u 2 Hx G +y G )^-mgu 3 = (39) 
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and the coefficient of degree one in u\, i.e. c5ml 

/ 2 \ dV 

2C (Cu 2 - mguiX G - mgu 5 y G ) qZ-x g u 3 u 2 

-AC {Cu 2 2 - mgu^xcj -^y G u 3 u 2 (40) 

/ 2 2 2 \ 3V 

-2 [2Cu 2 x G - 2mgu 4 x G + mgu 4 y G - mgu 5 x G y G j -^-mgx G u 2 

+ (Culx 2 G + 2Culy G - mgu A x G - mgu A x G y G ^ ^-mgu 3 = 

If V satisfies equations (|3"9*|) and Q4U[). then it is easy to prove that the 
determining equations Q32JI are identically satisfied by considering conditions 
(O, dSHD, ((SSI) as well. 

From (|39j) it is easy to obtain that V = ^(771,772,773) with 



Then, (|ifl|) becomes 



dV dV dV 

2mg(y G r} 2 - x G rj 3 )- h C— 773 - C—r] 2 = (42) 

d77i c% 0773 

Its characteristic curves are 

€i = Vi + ^fivcm + x G r) 2 ), £2 = V2 + »73 ( 43 ) 

Finally, we have that = ^(^1,^2) with \E' an arbitrary function of £i,£2> 
and consequently operator 

f = (44) 

is a generator of a Lie point symmetry for system (|25 |) -(|29 |) . Transforming 
(j43j) into the original unknown functions by using (|41|) , (j30|) , (|24j) , (|20|) yields 

6 " C I 2 



2p 2 + 2g 2 + r 2 ) + mg (x G a + y G /3)J 

x G a-y G /3\ 2 / x G /3 + y G aV 
rp z -q A -mg I +I2pg-m5 I 

6 = C 2 - 



m 2 g 2 (x G + y G ) 
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which correspond to the first integral of conservation of energy ((lip, and 
that derived by Kowalevski (|14|) . respectively 

Can other cases of integrability (viz integration by quadrature) be obtained 
by using our method? We leave the answer to a future paper. In [13], the 
application of our method led to an integrable case for a nonlinear system 
of three ordinary differential equations which does not possess the Painleve 
property. 

4 First integrals of the Kowalevski top 

We apply our method to the Kowalevski top itself which corresponds to the 
following conditions on the parameters: 

(1) A=B=2C 

(2) y G = z G = 0, x G > 

The condition on ug can be added without loss of generality 
(|lt)j) become 

p = rq/2 

q = —pr/2 + mgxGj /2C 

r = —rngxcP/C 

a = fir — ^q 

$ = jp — ar 

> 7 = aq- (3p 

The first integrals for the Kowalewski top are 

(1) conservation of energy, i.e. 

h = — [2p 2 + 2q 2 + r 2 J + mg x G a (46) 

(2) conservation of the vertical component of the angular momentum, i.e. 

h = C{2pa + 2q(3 + r 7 ) (47) 

(3) the length of the unit vertical vector, i.e. 

/ 3 = a 2 + /? 2 + 7 2 (=l) (48) 



. Then, system 



(45) 
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(4) the first integral derived by Kowalevski, i.e. 

/ 2 2 x G amg\ 2 ( x G (3mg\ 2 
h=\P ~q ^ — J + [2pq — 1 (49) 

If our method is applied to (|4*5|) . then all the first integrals can be obtained, 
apart from (J47|) which has all the unknown variables p, q, r, a, j3, 7 appearing 
in its expression. Let us observe that 

• j3 does not appear in I\ 

• 7 does not appear in both I\ and I4 

• p does not appear in I3 

In the following, we eliminate a,(3,j,p from system (|45|) one at a time. 

4.1 Eliminating a 

First we show a negative result: no first integral obtained. Let us assume 
that we do not know any of the first integrals. Therefore, we do not know 
a priori that none of the first integrals can be obtained by deriving a. We 
derive a from the fifth equation of system (|45j) . i.e. 

PI- $ 

a = 

r 

and obtain the following system of four equations of first order, and one of 
second order: 

"til = ( — 2CU1U3 — 2CM1U3U4 — 2CU1U2U4 + 3(7^2^3^5 

+2Cu2u\u§ — 2mgu\ii\XG + 2mgu\U4 : u^XG) /2Cus 
u 2 = (_Cu 3 ii4 + mgu 5 x G )/2C , . 

u 3 = (-mguxx G )/C 
u 4 = (u 2 u 3 )/2 

, U5 = (—U1U3U4 — il\U2 + U2U4Us)/u3 

with 

u\ = P, u 2 = q, u 3 = r, n 4 = p, u 5 = 7 (51) 

If we apply Lie group analysis to system (|5U|) , then we obtain a determining 
equation of parabolic type for V in two independent variables. Its charac- 
teristic curve is 

U5U3 + U1U2 
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which yields the following transformation 

s 5 - u 2 ui . s 5 -q/3 

U5 = viz 7 = (52) 

u 3 r 

with S5 a new unknown function of t. Then, system (|5U|) transforms into: 

iii = (SCuiulul - 2Cu\v^v\ - 2Cu\u\ - 2Cu\u 3 u\ 

—2CU1U2U3114 + 3CU2U3U5 + 2CM2W4U5 — 2mgu\u2U±XG 
—2mgu\U\u 3 x G + 2mgu\U4U^xa) /2Cu 3 
u 2 = {-Cului - mguiu 2 x G + mgu 5 x G )/2Cu 3 
u 3 = (-mguix G )/C 
U4 = u 2 u 3 /2 

U5 = (—2C111U2U4 — 3CuiU 3 U4 + 2CU2U4U5 + mgu\u2XQ 
-mgu 1 u 5 x G )/2Cu 3 

with 

u 5 = s 5 (54) 

If we apply Lie group analysis to system (|5Hj) , then we obtain a two dimen- 
sional Lie symmetry algebra generated by the following two operators: 

r, = §. (55) 

_, d d d d d 
T 2 = -t—+2ui~ hn 2 ^ \-u 3 - \-1i4- — 

OT OU\ OU2 OU 3 OU4 

+3n 57 |- (56) 
which in the original unknown functions correspond to: 

r ' = I < 57 > 

d d d d d 
ot op oq or oa 

This is a trivial finding. 
4.2 Eliminating [3 

We derive (5 from the third equation of system (|45|) . i.e. 
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and obtain the following system of four equations of first order, and one of 
second order: 



u x = (mgx G (uiU A - u 3 u 5 )/C 

u 2 = (-Cu 1 u 3 + mgx G u 5 )/2C 

u 3 = UiU 2 /2 

u 4 = {-CuiUi-mgx G U2U^)/mgx G 

M5 = {Cu\u 3 + mgx G U2Ui) /mgx G 



(59) 



with 



ui = r, u 2 = q, u 3 =p, 114 = a, u 5 = 7 (60) 

If we apply Lie group analysis to system (|59|) , then we obtain a determining 
equation of parabolic type for V in three independent variables. Its two 
characteristic curves yield the following transformations 



s 4 



Cu\ 



u 5 



2mg x G 

Cu 3 U\ + S6 

mg x G 



viz a 



viz 7 



s 4 - Cr 2 
2mg x G 

Cpr + s e 
mg x G 



(61) 



with S4 and sq new unknown functions of t. Then, system (|50j) transforms 
into: 

Hi = (— Cu\ — 2C111U3 + U1U4 — 2u 3 u^)/2C 

U2 = U5/2C 

u 3 = U1U2/2 (62) 

u 4 = -2u2(Cuiu 3 + u 5 ) 

U5 = u 2 {-2Cu\ + ua)/2 



with 

u 4 = S4, «5 = «6 (63) 

If we apply Lie group analysis to system (|62|) , then we obtain two first order 
partial differential equations for V: 



dV 

du 3 



dV 
8114 

dV 











(64) 
(65) 



dV 

du2 * du^ 

with V = V(u2, U3, U4). From it is easy to obtain that V = V{r),u 2 ) 
with 

7/ = 2Cu\ + u 4 (66) 
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Then, (|63)) becomes 

dV . „ dV 



ACu 2 — = (67) 
OU2 or] 



Its characteristic curve is 



£1 = 2CU 2 . + 77 (68) 



Finally, we have that V = ^(£1) with ip an arbitrary function of £, and 
consequently operator 

Ti = V(6)9 4 (69) 

is a generator of a Lie point symmetry for system 1)621) . Transforming (|68[) 
into the original unknown functions by using (|66|) . I|63|). (|61|). (|60|) yields 

£1 = - (2p 2 + 2g 2 + r 2 ) + m 5 x G a 

which is exactly the first integral of conservation of energy (|46[). In addi- 
tion, we have algorithmically derived that ()69() is a generator of a Lie point 
symmetry for system (|45|) . 

4.3 Eliminating 7 

We derive 7 from the second equation of system (|45[). i.e. 

C(2g + pr) 



7 



mgx G 



and obtain the following system of four equations of first order, and one of 
second order: 

111 = Cu 2 + 2mgu±XG)/AC 

112 = U1U3/2 

u 3 = -mgu 5 x G /C (70) 

U4 = (—2Cuiui — CU1U2113 + mgu^u^xo) /mgxG 

t ^5 = {2Cu\u 2 + Cu|u 3 - mgu 3 U4X G )/mgxG 

with 

it! = q, u 2 = p, u 3 = r, U4 = a, us = (3 (71) 

If we apply Lie group analysis to system l|7flj) , then we obtain a determining 
equation of parabolic type for V in three independent variables. Its two 
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characteristic curves yield the following transformations 

s 4 — Cu\ s 4 — Cq 2 

U4 = viz a - 



mgx G mgx G 

2Cu x u 2 + s 5 . 2Cpq + s 5 

«5 = viz p = (72) 

mgx G mgx G 

with s 4 and S5 new unknown functions of t. Then, system (|7U|) transforms 
into: 

ui = [ui(-2Cuf-Cul + 2u 4 )]/4C 

U2 = UiU 3 /2 

u 3 = (-2Cuiu 2 -u 5 )/C (73) 
u 4 = u 3 (Cuiu 2 + us) 

with 

U4 = s 4 , u 5 = s 5 (74) 

If we apply Lie group analysis to system ((73*)) , then we obtain two first order 
partial differential equations for V: 

dV dV dV 

u 3 - 4u 2 t^— + ±Cu 2 u 3 — = (75) 

OU 2 OU 3 OU4 

~ 9V „ „ dV „ /_ „ ,\ 9V 

8C-U2«5t; Cu 3 ii 5 — - + Gu 3 u 4 - Cu 2 ) — - = (76) 

ou 3 0U4 V / Otis 

with V = V(u2,u 3 ,U4,u^). From ((75]) it is easy to obtain that V = 
V(rji,V2,u 5 ) with 

rfr = Au 2 . + uf, 772 = Cu^ - tt 4 (77) 

Then, ((7""|) becomes 

2C« 3 u 5 75— + Cu 3 u§ Cu 3 ( u 4 - Cu\ ) — = 78 

OT]i OT]2 V ' C/M5 

Its characteristic curves are 

£1 = Ct?! - 2r, 2 , e 2 = ??I + «i (79) 

Finally, we have that V = ^(^i, £2) with \& an arbitrary function of £1, £ 2j 
and consequently operator 

ri = *(6,6)^ (so) 
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is a generator of a Lie point symmetry for system (|73|1. Transforming (|79jl 
into the original unknown functions by using (|77|). ifTljl . ((72J), (f7T|) yields 

£1 = - (2p 2 + 2g 2 + r 2 ) + mg x G a 

2 _2 x G amg\ 2 : x G f3mg\ 2 



which are exactly the first integral of conservation of energy (|46fl . and that 
derived by Kowalevski ()49j) . respectively. In addition, we have algorithmi- 
cally derived that (j80[) is a generator of a Lie point symmetry for system 

4.4 Eliminating p 

We derive p from the second equation of system (|45jl . i.e. 

mg'jXG - 2Cq 



P 



Cr 



and obtain the following system of four equations of first order, and one of 
second order: 

u\ = u\{—Cu 2 + 2mgu4 : XG)/^C 
u 2 = (CU1U3U4 + 2Ctt 1 ti 5 — mgu2U^XG)/Cu3 
u 3 = -mgu 5 x G /C (81) 
u 4 = -u\U2 + u s u 5 

U5 = {—2Cu\U2 — C113U4 + mgu\xQ)jCu 3 
with 

ui = q, u 2 = 7, u 3 = r, u 4 = a, u 5 = (3 (82) 

If we apply Lie group analysis to system (|81|) , then we obtain a determining 
equation of parabolic type for V in three independent variables. Its two 
characteristic curves yield the following transformations 

u 2 = y/s6 cos(2s5 — 2u 3 ) viz 7 = ■ v /sgcos(2s5 — 2r) 

= ^/sq sin(2s5 — 2^3) viz (3 = sin(2ss — 2r) (83) 
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with sq e S5 new unknown functions oft. Then, system ()81j) transforms into: 

iti(— Cu\ + 2mguiXG)/^C 
2m4V^cos(2u3 — 2tt 5 ) tan(2n 3 — 2u^)u 3 + U\ 
-V/S2 sin(2u3 — 2u^)mgxG / C 



ui 

U3 
U4 

U 5 



U2 cos(2-u 3 — 2tt 5 ) tan(2n3 — 2u^,)u 3 + u\ 
U2 cos(2n3 — 2u 5 ) — (C113U4 — mgu-zXQ) — 2Cu\U2 
+{Cu\U4 : + 2mgu2XG) tan(2«3 — 2115)7x3 /2Cu2U 3 



(84) 



with 



u 2 = s 6 , n 5 = s 5 (85) 

If we apply Lie group analysis to system (|84|) , then we obtain a determining 
equation of parabolic type for V in two independent variables. Its charac- 
teristic curve yields the following transformation 



u 5 



SS5 — U\ 

u 3 



VIZ 



7 : 



Sq cos (2 
si sin(2 



ss 5 



r 

ss 5 - q 



(86) 



with SS5 a new unknown function of t. Then, system (|84[) transforms into 



U2 
"5 



2ti4-^/u2 cos ((2ui + 2u| — 2u^)/u 3 ^ v,\ 

+ tan f(2«i + 2-Ug — 2u^)/u^ju 3 
\fu2~ sin ^(2iii + 2^3 — 2u§) / uzjmgxG / C 
— \fu2~ cos n2«i + 2u§ — 2u 5 )/u 3 

\fu2~ cos ( (2«| — 2-u 5 + 2«i) /ii 3 



«1 



+ tan ^(2«i + 2tig — 2-u 5 )/ti 3 )u3 



(87) 



[2u\ + 2u 5 



2ui)mgu2XG + Cuiu^u^ tan ^(2n| — 2-% + 2ui)/u 3 ^ 



+ {mgu 2 x G ~ cului)U3 



/2ctt 2 U3 



with 



ii 5 = ss 5 

If we apply Lie group analysis to system (|H7)l , then we obtain one first order 
partial differential equation for V: 



dV 

&U4 



2U4 



dV 
dm 







(89) 
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with V = V(u2, M4). Its characteristic curve is 

= u 2 + U4 (9°) 
Finally, we have that V = ip(£i) with ^ an arbitrary function of £1, and 
consequently operator 

ri = V(&)$ (9i) 

is a generator of a Lie point symmetry for system (|87|) . Transforming (|9Ujl 
into the original unknown functions by using ([8SJ). (|5B|l. (jHSJ), (|8"2*|) 

yields 

6 = q 2 + /? 2 + 7 2 

which is exactly the first integral of the length of the unit vertical vector 
(|48|) . In addition, we have algorithmically derived that (|91(l is a generator 
of a Lie point symmetry for system (|45|) 



5 Kepler problem 

In |33| . Nucci's method [HOI was used to find symmetries additional to those 
reported by Krause m his study of the complete symmetry group of the 
Kepler problem. A consequence of the application of Nucci's method was 
the demonstration of the group theoretical relationship between the simple 
harmonic oscillator and the Kepler problem. In |33| . polar coordinates were 
used, and Nucci's method was not applied to the three-dimensional case 
with the purpose of finding first integrals. We do it here by considering 
cartesian coordinates. 

The equations of motion of the Kepler problem are given by the following 
well-known three equations of second order 

x\ = — Lixi/( (x\ + x\ + X3) 3 / 2 ) 
< x 2 = -^x 2 /\Axl + xl + xlf/ 2 \ (92) 
x 3 = -fxx 3 /((xl + x\ + a: 2 ) 3 / 2 ) 

The first integrals for the Kepler problem are: conservation of energy E, 
conservation of angular momentum K, the Laplace-Runge-Lenz vector L. 
None of the unknowns x\, x 2 , X3, x\, x 2 , X3 are missing in the expression of 
E and the components of L. This is not true for the three components of 
K, i.e. 

K\ = x 3 x 2 - X3X2 (93) 
K 2 = x 3 xi - x 3 xi (94) 
K 3 = x\x 2 - x\x 2 (95) 
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Therefore, we can only obtain the three components of K using our method. 
However, neither E nor L are needed to reduce system ()92j) to a linear 
oscillator, as we show in the following. Let us transform system (j92|) into a 
system of six equations of first order 

W\ = w& 

W 2 = W 5 
W 3 = W 6 

< W4 = -^ Wl /(Awl + wl + wlf/ 2 \ (96) 

W 5 = -fJ 1 W 2 /((wl + W2+W%) 3 / 2 ) 

w 6 = -fiw 3 /(^(wl + wl + wl) 3 / 2 j 

with 

Wl = Xl, W 2 = X 2 , W 3 = X 3 , W4 = ±1, 1^5 = X 2 , Wq = X 3 (97) 

Consequently, the components of the angular momentum become 

K\ = w 3 w§ — wqw 2 (98) 
K 2 = w 3 W4 — w\Wq (99) 
K 3 = wiw§ — waw 2 (100) 

We choose one of the dependent variables to be the new independent variable 
y in order to reduce the order of system (|U6^1 by one jSOj. We take w 3 = y. 
Then, system (|96|) becomes the following non-autonomous system of five 
first order ordinary differential equations 

w[ = W4/WQ 

w' 2 = w 5 /w 6 
< w' A = -iiwi/{wQ{w 2 + wl + y 2 Y /2 
w' 5 = -nw 2 /(w 6 (wl + wj + y 2 ) 3 ^ 2 
k w 'e = ~lJ.y/(w & {wl + wl + y 2 f' 2 j 
with ' denoting differentiation with respect to y. Let us observe that 

• UI4 does not appear in K\ 

• W5 does not appear in K 2 

• wq does not appear in K 3 

We should remark that other variables are missing too. For example, w\ is 
also missing in K\. However, our method will yield the result whatever the 
choice of a missing variable. In the following, we eliminate w^,w^,wq from 
system one at a time. 
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5.1 Eliminating W4 

We derive from the first equation of svstem (|l(Jl|) . i.e. 

and obtain the following non-autonomous system of three equations of first 
order, and one of second order: 

' u'i 



u 2 



(i(u' 4 y - u A )/(ul{u 2 2 + u\ + y 2 fl 2 
-fiy/(u 3 (u 2 + uj + y 2 ) 3 / 2 ) 
ui/u 3 

-fiu 2 /{u 3 (u 2 2 + ul + y 2 f/ 2 ) 



with 



(102) 



(103) 



1i 4 = Wl, U 2 = W 2 , U 3 = W 6 , Ui = W 5 

If we apply Lie group analysis to system (|102|) , then after several reductions 
we obtain one first order partial differential equations for G 3 : 



8G 3 dG 3 
ui— h u 3 - 



du 2 dy ^ ^ 

with G 3 = G 3 (u\,u 2 ,u 3 ,y). Its solution is G3 = V'(Ci) with ifj an arbitrary 
function of 

6 = U3U2 - yui (105) 

Transforming (|1U5() into the original unknown functions by using (|97|) . H103|) 

yields 

£1 = X3X2 ~ x 2 x 3 

which is exactly the first component of the angular momentum (|93|) . 
5.2 Eliminating 105 

We derive w§ from the second equation of system (|101j) . i.e. 

w 5 = w' 2 w 6 

and obtain the following non-autonomous system of three equations of first 
order, and one of second order: 

' u'i = v{-UA + u' A y)/(u 2 {ul + u 2 + y 2 fl 2 ^ 

-mi/(Mu 2 + u 2 + y 2 fl 2 ) 

-^/({ui + ui + y 2 )*' 2 ) 
u 3 /u 2 
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with 

U 4 = W 2 , u 2 = We, U 3 = W4, U\ = W\ (107) 

If we apply Lie group analysis to system ()106|) , then after several reductions 
we obtain one first order partial differential equation for G 2 : 

dG 2 dG 2 , x 

u 3ir ^ + u 2 —^ = 108 
ou\ oy 

with G 2 = G 2 (u\,u 2 ,u 3 ,y). Its solution is G 2 = (f>(£ 2 ) with 4> an arbitrary 
function of 

6 = IJU3 ~ uiu 2 (109) 

Transforming (j!09[) into the original unknown functions by using (j97j) . (|107|) 
yields 

which is exactly the second component of the angular momentum ()94|) . 
5.3 Eliminating wq 

We derive wq from the first equation of svstem ljlOlj) . i.e. 

w 6 = — 
w\ 

and obtain the following non-autonomous system of three equations of first 
order, and one of second order: 

' u'i = ^u'i{-u 4 + u' A y)/(ul{ul+ul + y 2 f/ 2 ) 

u' 3 = -nmu'J^miul + ul + yy/ 2 ) , no . 

u' 2 = u 3 u'Ju\ 

u[ = -nu4u'J[ux{u\ + u\ + y 2 ) 3/2 ) 
with 

u 4 = wi, u 2 = w 2 , U 3 = w 5 , U\ = w 4 (111) 

If we apply Lie group analysis to system ()110|) , then after several reductions 
we obtain one first order partial differential equation for G 3 : 

0G 3 dG 3 
u 3 - hui- — = (112) 

OU 2 OU\ 



24 



with G3 = G3(ui,U2,U3,U4) with G3 = G 3 (ui, u?, U3, y). Its solution is 
G3 = 99(^3) with ip an arbitrary function of 

£3 = «1«2 - U 4 U 3 (113) 

Transforming Q113|) into the original unknown functions by using (|9~7)l , (|111|) 

yields 

£3 = i 2 ari - ±ix 2 

which is exactly the third component of the angular momentum (|94l) . 
Now let us derive w§, and W2 from (|98|). (|9*9*|) and (|1U0|> . i.e. 

w 5 = (114) 

u>4 = (115) 

y 

u; 2 = (116) 

6 

with 6) 61 6 new unknown functions of y. Substituting (|114|) . (|115|) . 
(j!16j) into IjlOlJI , and deriving wq from the first equation yields the following 
system of three equations of first order, and one of second order: 

u'l = (fj,u 2 (-ul + 3ulu' 4 y -3U4U4V + u'ly % )^ / {[u\y 2 

-2u 3 uiU4,y + u\u\ + u\y 2 + n|n|) 3 / 2 ) 
u> 3 = " ' ( 117 ) 

«2 = 

ui = 
with 

■u 4 = wt, u 3 = 6, ^2 = 6> "1 = 6 

It is easy to show that system (j!17|l admits an eleven-dimensional Lie sym- 
metry algebra. In fact, the first equation of (|117|) itself admits a Lie symme- 
try algebra of dimension eight, which means that it is is linearizable through 
a point transformation Thus, we have reduced the equations of motion 
of the Kepler problem to the harmonic oscillator [33] , |35| by using Lie group 
analysis. 

6 Final comments 

We have found that Lie group analysis yields the first integrals admitted 
by any system of ordinary differential equations if the method developed by 
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Nucci |3Uj is applied, the only limitation being the absence of at least one 
of the unknowns in each first integral. 

Is it possible to obtain all of the first integrals by means of Lie group analysis? 
Also, what is the link between Painleve method and Lie group analysis |34j ? 
In addition, can Lax pairs be found by Lie group analysis? So far these are 
open questions that we hope to address in future work. 
Let us conclude by underlining that the application of Nucci's method to the 
Kowalevski top have led us to understand how first integrals can be found 
by using Lie group analysis. In 1984, Cooke [I] wrote "Kowalevskaya's 
work is an ingenious application of mathematics to a system of equations 
of great mathematical interest . . . but since the case to which it applies is 
rather special, the details of her arguments are no longer worth troubling 
about". About the same time, a revival of interest into integrable problems 
of mechanics has led to numerous papers on the Kowalevski top. Just to 
cite a few, see [HI], EE] IE] P |S| [E| and the entire no. 11 issue of the 
Journal of Physics A 34, (2001). 
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